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Abstract 

Background: Hox genes encode master regulators of regional fate specification during early metazoan development. 
Much is known about the initiation and regulation of Hox gene expression in Drosophila and vertebrates, but less is 
known in the non-arthropod invertebrate model system, C. elegans. The C. elegans Hox gene lin-39 is required for 
correct fate specification in the midbody region, including the Vulval Precursor Cells (VPCs). To better understand lin-39 
regulation and function, we aimed to identify transcription factors necessary for lin-39 expression in the VPCs, and in 
particular sought factors that initiate lin-39 expression in the embryo. 

Results: We used the yeast one-hybrid (Y1 H) method to screen for factors that bound to 1 3 fragments from the lin-39 
region: twelve fragments contained sequences conserved between C. elegans and two other nematode species, while 
one fragment was known to drive reporter gene expression in the early embryo in cells that generate the VPCs. Sixteen 
transcription factors that bind to eight lin-39 genomic fragments were identified in yeast, and we characterized several 
factors by verifying their physical interactions in vitro, and showing that reduction of their function leads to alterations 
in lin-39 levels and lin-39::GFP reporter expression in vivo. Three factors, the orphan nuclear hormone receptor NHR-43, 
the hypodermal fate regulator LIN-26, and the GATA factor ELT-6 positively regulate lin-39 expression in the embryonic 
precursors to the VPCs. In particular, ELT-6 interacts with an enhancer that drives GFP expression in the early embryo, 
and the ELT-6 site we identified is necessary for proper embryonic expression. These three factors, along with the 
factors ZTF-17, BED-3 and TBX-9, also positively regulate lin-39 expression in the larval VPCs. 

Conclusions: These results significantly expand the number of factors known to directly bind and regulate lin-39 
expression, identify the first factors required for lin-39 expression in the embryo, and hint at a positive feedback 
mechanism involving GATA factors that maintains lin-39 expression in the vulval lineage. This work indicates that, as in 
other organisms, the regulation of Hox gene expression in C. elegans is complicated, redundant and robust. 
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Background 

Hox genes encode evolutionarily conserved homeodomain- 
containing transcription factors that pattern cells and 
tissues along the anterior-posterior body axis during 
metazoan development (for review, see [1,2]). Hox pro- 
teins perform this function by serving as master regulators 
of expression of batteries of genes that impart identity to a 
cell [3,4], and the precise regulation of Hox protein activ- 
ity is therefore vital for proper development. Due to their 
central and conserved role in regional identity and fate 
specification during metazoan development, the mecha- 
nisms underlying the initiation, maintenance and modula- 
tion of Hox gene expression have been intensively studied 
[5-8]. In Drosophila, an elaborate network of maternally 
supplied and zygotically expressed transcription factors 
act to initiate Hox gene expression properly in the syncyt- 
ial early embryo, while in the cellularized vertebrate em- 
bryo, secreted signaling molecules and growth factors are 
employed to coordinate precise Hox gene expression in 
time and space [9,10]. In both vertebrates and Drosophila, 
once Hox gene expression is initiated, it is maintained by 
the Trithorax and Polycomb groups of chromatin regula- 
tory proteins [11,12], and can be further modulated by 
extracellular signaling pathways, Hox protein autoregula- 
tion and cross-regulation, and other mechanisms [5,13]. 

As in other metazoans, Hox genes are essential during 
the development of the nematode Caenorhabditis elegans 
[14]. C. elegans has only six Hox genes, present in a dis- 
persed cluster [15,16]. Three Hox genes, ceh-13, nob-1, 
and php-3, are required for proper embryonic develop- 
ment [17-19]. The other three Hox genes, lin-39, mab-5, 
and egl-5, appear to be required only during post- 
embryonic development, however their expression begins 
in the embryo [20-22] . Very little is known about the initi- 
ation of Hox gene regulation in C. elegans in the embryo 
[23], although subsequent regulation of Hox gene expres- 
sion in larval development by Polycomb and Trithorax- 
related proteins, signaling pathways, other transcription 
factors, microRNAs, and Hox proteins themselves, have 
all been noted [12,23-40]. Elucidating the mechanisms by 
which Hox gene expression is initiated and regulated in 
nematodes will broaden our understanding of this import- 
ant class of developmental regulators across a larger range 
of animal phyla, giving us further insight into their use 
during the evolution of animal diversity and their function 
in gene regulatory networks controlling pattern formation. 

Our laboratory and others have studied the function of 
the Hox gene lin-39 during nematode larval develop- 
ment, in particular during formation of the vulva, which 
is part of the hermaphrodite egg-laying apparatus. Vulval 
development begins in the first larval stage (LI) when 
the twelve ventral hypodermal blast cells, P1-P12 (P 
cells), divide to generate posterior daughters known as 
Pn.p cells [41]: the central six Pn.p cells, P3.p-P8.p, 



become Vulval Precursor Cells (VPCs) [42]. During the 
third larval stage (L3) the action of Wnt, Ras and Notch 
extracellular signaling pathways induces the VPCs to 
adopt distinct cell fates in the pattern 3°-3°-2°-l 0 -2 0 -3°, 
where the cells P5.p - P7.p adopt 1° and 2° (vulval) fates 
and divide to generate 22 cells that form the vulval 
opening, while P3.p, P4.p and P8.p adopt the non-vulval 
3° fate, which is to divide once and fuse with the sur- 
rounding syncytial hypodermis (reviewed in [43,44]). 
The Hox gene lin-39 encodes a Deformed/Sex combs 
reduced ortholog expressed in the midbody region, in- 
cluding the six VPCs [20,22]. lin-39 acts twice during 
vulval development, lin-39 is first required to generate 
the VPCs; in lin-39 null mutants, the VPCs fuse with the 
hypodermis during the LI stage, causing a Vulvaless 
(Vul) phenotype [20,22]. Little is known about the regu- 
lation of lin-39 expression at this time in development. 
lin-39 is also required at the time of VPC fate specifica- 
tion in the L3; loss of lin-39 activity at this time leads to 
defects in VPC fate specification [28,45]. At this later 
time, LIN-39 acts downstream of RTK/Ras and Wnt 
extracellular signaling pathways [26,28,46,47]. 

Trans-acting factors regulating lin-39 expression have 
been identified previously using a variety of methods 
including forward and reverse genetic analysis, evolu- 
tionary conservation, and transgenic reporter analysis. 
Trans-acting factors regulating lin-39 expression during 
vulval development include the RTK/Ras pathway tran- 
scriptional effectors LIN-1 and LIN-31 [28,46-48], the 
Gli family member TRA-1 which acts downstream from 
the sex determination pathway [37], the zinc finger pro- 
tein SEM-4 [49], the novel protein LIN-25 [47], several 
chromatin regulators [25,32,34,48,50], and LIN-39 itself 
[28,47]. Direct binding to sites within the lin-39 genomic 
region has been established in the case of LIN-1 (and 
LET-418, with which it interacts), LIN-31, LIN-39 and 
TRA-1 [37,46-48]. 

We previously undertook to identify cis-acting sites 
regulating lin-39 expression [47]. Due to the large size of 
the lin-39 genomic region (~28 kb) we inserted fragments 
of lin-39 genomic DNA upstream of an enhancerless GFP 
reporter. By that method, we identified three functional el- 
ements: a 340 bp upstream fragment (JW3.9) that directs 
expression in P cells in the embryo and in their larval de- 
scendants including the VPCs, a 247 bp site from the first 
lin-39 intron that directs expression in larval ventral cord 
neurons, and a 1.3 kb promoter fragment (JW5) that 
drives expression in P6.p at the time of vulval induction. 
Expression from the last element is dependent on Ras 
pathway function and we showed that LIN-1, LIN-31 and 
LIN-39 directly bind this cis-regulatory module. Using an 
alternative approach, Kuntz et al. used phylogenetic ana- 
lysis to identify a number evolutionarily-conserved regions 
in the lin-39 genomic region, and showed that several of 
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these sites also drove expression in certain cells or tissue 
when inserted upstream of an enhancerless GFP reporter 
[51]. These approaches both require that the site in ques- 
tion be able to mediate proper transcription activation of 
the reporter gene in vivo, and this requirement may lead 
to false negative results. To circumvent this issue, here we 
use the yeast one-hybrid assay, in which transcription fac- 
tors that directly bind to a site of interest are identified, 
and then the function of these factors in gene regulation 
in vivo can be assayed [52] . 

In our previous analysis of the lin-39 genomic region, 
we identified a number of short DNA sequences in the 
lin-39 region that were strongly conserved between C. 
elegans and two other nematode species ([47], Supple- 
mental material). Starting with these evolutionarily con- 
served elements and the 340 bp JW3.9 fragment, we 
used manual [53] and robotically-assisted "enhanced 
Y1H" (eYIH) [54] screens to identify 16 C. elegans tran- 
scription factors that bind to specific lin-39 genomic 
DNA fragments. Seven of these factors were character- 
ized further to determine a role in lin-39 regulation dur- 
ing vulval development. We found that the orphan 
nuclear receptor NHR-43, the hypodermal fate regulator 
LIN-26 and the GATA factor ELT-6 positively regulate 
lin-39 expression in the embryo and may play a role in 
initiation of lin-39 in the vulval lineage. In the larva, 
NHR-43, LIN-26, ELT-6, the zinc finger proteins ZTF-17 
and BED-3, and the T box factor TBX-9 positively regu- 
late lin-39 expression in the VPCs. Interestingly, we 
previously showed that the adjacent GATA factor genes 
egl-18 and elt-6 are downstream targets of LIN-39 in the 
larva VPCs [55]. Combined with our current result that 
ELT-6 binds to and regulates lin-39 expression in the 
embryo, this suggests that EGL-18/ELT-6 and LIN-39 
may form a positive feedback loop to initiate and main- 
tain lin-39 gene expression during embryonic and larval 
life to ensure proper VPC fate specification. 

Results 

Identification of transcription factors that bind to lin-39 
genomic regions using yeast one-hybrid screens 

We previously identified cis and trans-acting factors that 
regulate lin-39 expression in response to RTK/Ras sig- 
naling in the larval VPCs [47]. To further understand 
Hox gene lin-39 regulation in C. elegans, we sought to 
identify transcription factors that bind to and regulate 
lin-39 expression, with an emphasis on 1) regulation in 
the vulval precursor cells P3.p - P8.p (VPCs) and their 
descendants in the larva, and 2) expression in the pre- 
cursors to the VPCs, the P cells P3 - P8, in the embryo. 
We used the yeast one-hybrid approach, in which DNA 
sequences from the gene of interest are used as 'bait' to 
screen for factors that can bind these sequences in yeast 
leading to activation of reporter gene expression [53]. 



Previously, in addition to identifying large genomic re- 
gions that drove GFP reporter expression in lin-39- 
expressing cells, we also identified 31 short DNA elements 
(<50 bp) located upstream, downstream and in introns 
that are conserved in lin-39 from the species C. elegans, C. 
briggsae and C. remanei [47] (Additional file 1: Figure SI). 
Reasoning that some of these evolutionarily conserved re- 
gions (ECRs) may be binding sites for transcription factors 
that regulate lin-39 expression, we used twelve small lin- 
39 genomic regions (range 150-460 bp) that each encom- 
passed one or more ECRs in yeast one-hybrid screens 
(Table 1, Additional file 1: Figure SI). We also used the 
element JW3.9, a 340 bp fragment found 7.4 kb upstream 
of lin-39 that drives GFP expression in the P cells of the 
embryo [47] , as this fragment may bind factors responsible 
for initiation of lin-39 expression (Additional file 1: Figure 
SI). These thirteen DNA fragments, which together repre- 
sent ~ 13% of the genomic region between lin-39 and its 
neighboring genes, were used as separate 'baits' (Table 1). 

We used two yeast one-hybrid assay procedures. One 
screen was performed by transformation of the thirteen 
bait strains with a library of 755 plasmids that each ex- 
press one C. elegans transcription factor fused to the ac- 
tivation domain from the yeast transcription factor 
GAL4 [53,56]. Two other screens were performed using 
a robotically-assisted mating assay in which the thirteen 
haploid bait strains were mated to a collection of 936 
strains, each of which expresses a single C. elegans tran- 
scription factor fused to the GAL4 activation domain 
[54]. Factors identified as positive from the primary 
screens were retransformed manually and tested a 
second time; only factors that showed a reproducible 



Table 1 Seven transcription factors interact with lin-39 
genomic DNA in yeast one-hybrid assays 


Fragment 


ECRs 


Size 


TFs bound 


YF1 


ECR1 


342 


NHR-43 


YF2 


ECR2 


311 


ALR-1 


YF3 


ECR4 


253 




YF4 


ECR7 - 10 


372 


ZTF-17, LIN-26 


YF5 


ECR11,12 


298 




YF6 


ECR16 


304 




YF7 


ECR17 


158 




YF8 


ECR18 - 20 


331 




YF9 


ECR21 - 23 


257 


TBX-9 


YF10 


ECR24 - 26 


319 


BED-3 


YF11 


ECR27,28 


254 




YF12 


ECR29 - 33 


455 




JW3.9 




338 


ELT-6 



The table indicates the 13 fragments (YF1-12, JW3.9) used in yeast one-hybrid 
screens, the evo I utionarily-con served regions (ECRs) each contains, the size of 
the fragment used, and the transcription factors identified as binding to each 
fragment. 
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interaction in yeast were considered true positives. In 
total, 18 interactions were found in yeast between 16 tran- 
scription factors and eight different fragments (site bound 
by each factor is shown in parenthesis): NHR-43 (YF1), 
ALR-1 (YF2), and ZTF-17 (YF4) were identified through 
the library transformation screen, while 13 factors were 
identified via the robotically-assisted mating screens: 
ODR-7, TBX-39 and TBX-40 (YF1); LIN-26 (YF4); TBX- 
11, TBX-39 and EGL-43 (YF7); FLH-1 and NHR-111 
(YF8); TBX-9 and the protein encoded by B0238.ll 
(YF9); BED-3 and FLH-1 (YF10); DMD-3 and ELT-6 
(JW3.9) (Additional file 2: Figure S2 and Additional 
file 3: Figure S3, Additional file 4: Table S2 and S3, 
summary in Additional file 4: Table S4). 

We describe here our characterization of seven tran- 
scription regulators (Table 1) chosen because these fac- 
tors had either known expression in Ziw-39-expressing 
cells or a phenotype affecting a Zi'w-39-regulated process 
[57], or because our preliminary data showed an effect 
on lin-39 reporter expression in vivo. For each factor we 
carried out the following analyses. To validate the yeast 
interactions we expressed and purified each factor from 
bacteria and assayed binding to the appropriate sites 
in vitro (Figure 1, Additional file 5: Figure S6). To deter- 
mine if these factors regulate lin-39 expression in vivo in 
the vulval precursor cells during larval life, we reduced 
function for each transcription factor and examined lin- 
39 expression using an integrated transcriptional lin-39:: 
GFP reporter (dels4) that contains 250 kb of genomic 
DNA around the lin-39 locus [46] (Figures 2 and 3). We 
also examined lin-39 expression in RNAi-treated and 
mutant animals at the L3 stage by qRT-PCR (Figure 4). 
The strain containing the integrated lin-39::GFP reporter 
was also used to examine the effects of reduced tran- 
scription factor activity on lin-39 expression in the em- 
bryonic P cells P5 - P8, which divide in the larva to 
generate /««-39-expressing vulval precursor cells and 
ventral cord neurons [46] (Tables 2 and 3, Figure 5). Fi- 
nally, we examined the effect of reduction of factor func- 
tion in a //«-39-sensitized background on the fusion of 
the vulval precursor cells at the L2 stage (Table 4). 
Below we describe our results that show that six of the 
factors identified in the yeast one-hybrid screens regu- 
late lin-39 expression in the larval vulval precursor cells, 
while three of them also regulate lin-39 expression in 
the embryo. 

Orphan nuclear hormone receptor NHR-43 positively 
regulates lin-39 expression in the embryo and larva 

The nhr-43 gene encodes an orphan nuclear hormone 
receptor, and nhr-43::GFP reporter expression is present 
from late embryo to adult in hypodermis, excretory cells, 
posterior intestine cells, and two head neurons [58,59]. 
In yeast, NHR-43 bound to YF1, a 342 bp fragment 



located 8.6 kb upstream of lin-39 (Additional file 1: 
Figure SI and Additional file 2: Figure S2), and this 
binding was verified in vitro using NHR-43 protein puri- 
fied from E. coli (Figure 1A). Some nuclear hormone re- 
ceptors are known to bind sites containing the sequence 
TGAC [60], and there is a TGAC site in YF1. While the 
wild type YF1 sequence competed for NHR-43 binding 
in vitro, the same YF1 fragment with the TGAC site mu- 
tated no longer competed effectively (Figure 1A), sug- 
gesting that NHR-43 may bind to YF1 through this 
putative NHR binding site. 

When RNAi was performed on worms carrying the 
lin-39::GFP transcriptional reporter to reduce nhr-43 
function, their progeny showed a decrease in the num- 
ber of animals expressing wild type levels of GFP in the 
vulval precursor cells at the L3 stage (Figures 2B and 
3A). qRT-PCR on nhr-43(tml381) mutant animals [57] 
at the L3 stage also showed a decrease in lin-39 expres- 
sion compared to control animals (Figure 4A). We also 
found that fewer embryos derived from hermaphrodites 
treated with nhr-43 RNAi showed lin-39::GFP expression 
in P5 - 8 compared to control animals (Figure 5B, 
Table 2), indicating that nhr-43 positively regulates lin- 
39::GFP in the embryonic P cells. This embryonic defect 
could explain the reduction in GFP expression in the 
VPCs in the larva, however a decrease in lin-39::GFP ex- 
pression was also observed if nhr-43 RNAi was per- 
formed on LI worms and then the same animals were 
observed in the L3 stage (51% of nhr-43(Ll RNAi) ani- 
mals showed wild type expression in all VPCs versus 
76% of RNAi control, p < 0.001). In summary, orphan 
nuclear hormone receptor NHR-43 binds to a site lo- 
cated far upstream of the lin-39 initiation codon, and 
nhr-43 function is required for wild type levels of lin-39 
expression in the embryo and larva in cells that will par- 
ticipate in vulval development. 

The Aristaless homolog ALR-1 may regulate larval lin-39 
expression 

alr-1 encodes a factor homologous to the products of 
the Drosophila aristaless and mouse ARX (aristaless- 
related) genes, and antibody staining and GFP reporter 
analysis have shown alr-1 expression in multiple neu- 
rons and hypodermal cells from embryo to adult [61,62]. 
In yeast, ALR-1 bound to YF2, a 311 bp DNA fragment 
located 6.4 kb upstream of lin-39 containing the evolu- 
tionary conserved sequence ECR2 (Additional file 1: 
Figure SI and Additional file 2: Figure S2). Recently, the 
modENCODE project identified genomic binding sites 
for a number of C. elegans transcription factors, includ- 
ing ALR-1, using chromatin immunoprecipitation from 
larval animals [40,63]. Examination of this data shows 
binding of ALR-1 to multiple sites in the lin-39 promoter 
region, including the ECR2 site. Bacterially expressed 
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Figure 1 LIN-39 regulators identified in yeast bind to lin-39 promoter regions in vitro. A - F) Gel mobility shift assays with proteins expressed 
and purified from E. coll Arrowhead indicates free probe; arrow indicates protein DNA complexes. For each panel, the top line identifies the labeled 
probe used, the bottom line shows the amount of protein added in each lane; middle lines (panels A, B, D and F) indicate the identity and amount of 
competitor added. A) NHR-43 binds to YF1 (342 bp). Lanes 3-5 contain cold wild type YF1 as competitor, while lanes 6-8 contain cold YF1 with the 
TGAC site mutated as competitor; B) ALR-1 binds to ECR2 (40 bp); this binding is competed by wild type ECR2 but not by a scrambled oligonucleotide 
with the same nucleotide composition (ECR2S); C) ZTF-17 binds to YF4-2 (186 bp) and YF4-4 (1 10 bp); D) ZTF-17 binding to YF4-2 is competed by 
YF4-4 but not YF4-3; E) LIN-26 binds to YF4-3-1 (51 bp); F) LIN-26 binding to YF4-3-1 is competed by YF4-3-1 (51 bp) but not YF4-3-2 (52 b); 
G) Fragment YF4 with ECRs 7-1 0 is shown above, with smaller subfragments diagrammed below. Shading indicates fragments that were bound by 
ZTF-1 7 and/or LIN-26 in yeast one-hybrid assays and in vitro. 



ALR-1 protein bound to a 40 bp DNA fragment encom- 
passing ECR2 in vitro, and this binding was not com- 
peted by 100 fold excess of an oligonucleotide of the 
same base composition and length but scrambled se- 
quence (Figure IB). Therefore, ALR-1 binds a site in the 
upstream region of lin-39 in both in vitro and in vivo 
assays. 

We did not observe a significant effect of alr-1 RNAi 
on lin-39 or lin-39::GFP expression in either the embryo 



or larva in the cells that will give rise to the vulva 
(Figures 2C, 3B, and 4A, Table 2). However, we did ob- 
serve an effect of alr-1 RNAi on a lin-39 mediated 
process in the larval VPCs. During wild type develop- 
ment, the cells P3.p - P8.p are born in the LI stage and 
because they express lin-39 they do not fuse with the sur- 
rounding hyp7 syncytium as more anterior and posterior 
Pn.p cells do [20,22]. Later, in approximately 50% of ani- 
mals, P3.p can fuse with hyp7 in the L2 stage [42]. Thus 
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Figure 2 Seven transcription factors affect lin-39::GFP expression in the VPCs at early L3 stage. GFP expression in the VPCs P5.p - P8 p 
from smg-1; delS4 {lin-39::GFP) animals treated for RNAi of individual transcription factors identified in yeast (panels B - H), or given control RNAi 
(empty vector; panel A). Early L3 stage animals are shown; anterior is left, ventral is down. These animals also express ajm-1::GFP, which outlines 
cell junctions of hypodermal cells and the pharynx (bright staining in the anterior seen in most panels). RNAi for lin-26, bed-3, and tbx-9 causes 
severe embryonic lethal and larva arrest phenotypes, so RNAi for these genes was performed by feeding newly-hatched L1 larvae on RNAi bacterial 
lawns and examining GFP expression in these same animals at the L3 stage. For all other genes, RNAi treatment was carried out on PO animals, and 
their F1 progeny were examined as L3 animals. All pictures were taken under the same exposure. 



wild type animals have either five or six VPCs at the time 
of vulval induction in the late L2 stage. In lin-39 mutant 
animals, or in animals in which Wnt signaling is compro- 
mised, additional VPCs are seen to adopt this Fused fate 
[26,28,45] . To examine the role of potential lin-39 regula- 
tors in this process we used RNAi to reduce their func- 
tion in a sensitized strain containing the hypodermal 
junction marker ajm-l::GFP and the temperature sensi- 
tive mutation Un-39(n709ts) and examined the fusion of 
the VPCs at the L2 stage as before [55]. We found an in- 
creased number of larva with VPC fusion defects in lin- 
39(ts) alr-l(RNAi) compared to lin-39(ts) control animals 
(Table 4). Therefore, although we were unable to detect 
an effect of loss of ALR-1 function on lin-39 expression 



in the vulval cells, we did see a weak effect on a LIN-39 
dependent process, leaving open its role as an in vivo 
regulator of lin-39 expression. 

The zinc finger protein ZTF-17 positively regulates lin-39 
expression 

ztf-17 encodes a zinc finger transcription factor, and a 
ztf-17::GFP transcriptional reporter is expressed in head 
muscle, pharynx, and the ventral nerve cord [64]. In 
yeast, ZTF-17 bound to YF4, a 372 bp DNA fragment 
located 2 kb upstream of lin-39 that contains ECRs 7-10 
(Additional file 1: Figure SI and Additional file 2: Figure S2). 
To identify a smaller ZTF-17 binding region, we divided 
YF4 into smaller fragments, and showed by one-hybrid 
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Figure 3 Seven transcription factors affect lin-39::GFP expression in the VPCs at early L3 stage, smg-1; dels4(lin-39::GFP) animals were 
treated for RNAi of individual transcription factors (panels A - G, dark bar), or control RNAi (panels A - G, light bar) as in Figure 2. Early L3 stage 
animals were photographed at the same exposure, and pixel counts in each VPC were determined using ImageJ (>20 animals for each strain). 
Bars show mean GFP pixel count in each cell with standard deviation. '*' indicates P-value < 0.05. '**' indicates P-value < 0.001, compared to 
control animals. 



assays that ZTF-17 interacts with YF4-4, a 110 bp subfrag- 
ment that overlaps ECR10 (Figures 1C and E, Additional 
file 6: Figure S4). in vitro binding and competition assays 
with bacterially-expressed ZTF-17 protein validated this 
result, showing that ZTF-17 binds to fragment YF4-4, but 
not the adjacent fragment YF4-3 Figure 1C and D). 

ztf-17 RNAi performed on lin-39::GFP worms caused a 
decrease in the number of progeny animals showing wild 
type levels of lin-39::GFP expression in the VPCs, 



(Figures 2D and 3C), and qRT-PCR analysis of these ztf- 
17(RNAi) animals showed a statistically significant de- 
crease in lin-39 at the L3 stage (Figure 4A). lin-39::GFP 
expression in the VPC parent cells, P5 - P8, was not al- 
tered in embryos derived from ztf-17(RNAi) mothers 
(Table 2), suggesting that ZTF-17 may be a larval regula- 
tor of lin-39 expression. In summary, in vitro and yeast 
one-hybrid analyses indicate that the zinc finger protein 
ZTF-17 binds to a 110 bp fragment located upstream of 
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Figure 4 Regulation of //n-39 levels by transcription factors in vivo (qRT-PCR). A) Decrease in lin-39 transcript levels in nhr-43, lin-26, tbx-9, 
bed-3 mutant and ztf-17(RNAi) animals. B) Decrease in lin-39 transcript levels when activity of both elt-6 and egl-18 is reduced. C) Increase in lin-39 
transcript levels when either elt-6 or egl-18 was overexpressed from the heat shock promoter. All analyses were done on L3 stage animals. The 
mean values for each genotype were obtained from at least two independent experiments and normalized to the housekeeping gene, gpd-2, as 
internal standard. The data was analyzed with unpaired t-test compared to the appropriate control. '*' indicates P-value < 0.05. 



lin-39, and positively regulates lin-39 expression in the 
vulval precursor cells during larval life. 

LIN-26, which is required for maintenance of hypodermal 
cell fates, may positively regulate lin-39 expression in the 
embryo and larva 

lin-26 encodes a zinc-finger protein which is expressed 
in all hypodermal cells based on antibody staining 
[65,66]. In lin-26 mutants, cell fate transformations from 
hypodermal to neuronal fate occur in many cells, includ- 
ing the VPCs [65,66]. LIN-26 also bound to YF4 in the 
yeast one-hybrid assay (Additional file 1: Figure SI and 
Additional file 7: Figure S5), and further analysis showed 



an interaction with YF4-3, a 103 bp DNA fragment that is 
distinct from the fragment bound by ZTF-17 (Figure 1G, 
Additional file 5: Figure S5). We divided fragment YF4-3 
into two smaller fragments, and in vitro DNA binding and 
competition experiments showed that LIN-26 bound to 
YF4-3-1, a 51 bp DNA fragment that does not contain an 
evolutionarily conserved region (Figure IE and F). 

To assay lin-26 regulation of lin-39 expression in vivo, 
lin-26 RNAi was performed on newly hatched lin-39:: 
GFP larvae and GFP levels in VPCs P5.p-P8.p were ex- 
amined at the L3 stage in the same animals; this 'LI 
feeding' was performed because of the embryonic lethal 
phenorype caused by lin-26 RNAi. lin-26 RNAi caused a 
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Table 2 Decreases in lin-39::GFP expression in embryonic 
P cells in transcription factor RNAi animals 



RNAi 


N 




% eggs with WT 
expression in P5-P8 


Control 


46 




100% 


nhr-43 


40 




80%* 


alr-1 


36 




100% 


ztf-17 


36 




100% 


lin-26 


36 




86%* 


tbx-9 


37 




100% 


bed-3 


36 




94% 


elt-6 


51 




82%* 


egl-18 


37 




100% 


smg-1; lin-39::GFP animals were grown on bacteria expressing dsRNA for each 
transcription factor or the control vector and GFP expression in the P cells (P5 - 
P8) of embryos laid by these animals was examined. The percentage of animals 
having a wild-type intensity of GFP in all four P cells based on visual observation 
is indicated. '1ST = number of embryos scored. '*' indicates P < 0.05 compared to 
control (Fisher's exact test). 


weak but significant decrease in expression of lin-39:: 
GFP in P5.p and P6.p at the L3 stage (Figures 2E and 
3D), and qRT-PCR performed on lin-26(nl56) mutant 
animals [65] also showed a small but significant decrease 
in lin-39 at the early L3 stage (Figure 4A). Finally, we in- 
vestigated lin-26 regulation of lin-39 in the embryo, and 
found that fewer embryos expressed lin-39::GFP in the P 
cells from mothers treated with lin-26 RNAi, compared 
to control embryos (Figure 5C; Table 2). Consistent with 
effects on lin-39 expression, lin-26(RNAi) performed on 
newly hatched LI larvae in a sensitized lin-39 back- 
ground caused a defect in VPC fusion at the L2 (Table 4). 


Table 3 Decreases in pJW3.9::GFP expression in elt-6 and 
egl-18 mutants in embryonic P cells 


Background 


RNAi 


N 


% eggs with WT 
expression in P5-P8 


Wild-type 




50 


98% 


Wild-type 


Control 


36 


94% 


Wild-type 


elt-6 


50 


78%* 


Wild-type 


egl-18 


36 


91% 


elt-6(gk723) 


Control 


34 


56%* 


elt-6(gk723) 


egl-18 


36 


58%* 


egl-18(n162) 


Control 


38 


68%* 


egl-18(nl62) 


elt-6 


36 


69%* 


egl-18(ga97) 


Control 


20 


60%* 


egl-18(ga97) 


elt-6 


24 


54%* 



Animals containing pJW3.9::GFP in either wild type, e/f-6 mutant (gk723) or 
egl-18 mutant (n/62 and ga97) backgrounds were synchronized and given the 
indicated RNAi treatment. pJW3.9::GFP expression in P5 - P8 was scored at 
bean stage in the embryos laid by these animals. 'N' = number of embryos 
scored. The percentage of animals having a wild-type intensity of GFP in all 
four P cells based on visual observation is indicated in column four. '*' indicates 
P < 0.05 compared to control (Fisher's exact test). 



Thus, we have identified a binding site in the lin-39 up- 
stream region for the general hypodermal transcription 
factor LIN-26, and our results suggest that LIN-26 acts 
to positively regulate lin-39 expression in both the em- 
bryo and larva. 

The T-box protein TBX-9 binds in lin-39 intron 2 and 
positively regulates larval lin-39 expression in the VPCs 

tbx-9 encodes a T-box transcription factor, and a tbx-9:: 
GFP translational reporter is expressed in lateral and 
ventral hypodermal cells, intestine and muscle [67]. tbx- 
9 mutants have a disorganized body shape beginning in 
the embryo, which was attributed to defects in hypoder- 
mal cells and body wall muscles [67]. In yeast, TBX-9 
bound to YF9, a 257 bp DNA fragment located in the 
large second intron of lin-39 that contains ECRs 21-23 
(Additional file 1: Figure SI and Additional file 2: Figure S2), 
and this binding was recapitulated in vitro with purified 
TBX-9 protein (Additional file 5: Figure S6). We divided 
the YF9 fragment into four smaller fragments (A-D), and 
found that bacterially purified TBX-9 bound best to the 
79 bp subfragment C, and that this fragment could also 
compete for binding of TBX-9 to YF9 (Additional file 5: 
Figure S6). 

To test TBX-9 as a regulator of lin-39 in vivo, tbx-9 
RNAi was performed on LI lin-39::GFP worms, and we 
observed a decrease in the number of animals with wild 
type levels of GFP expression in P5.p-P8.p in the same 
animals at the L3 stage (Figures 2F and 3E). qRT-PCR 
analysis on tbx-9(ok2473) mutant animals [57] also 
showed a decrease in lin-39 levels compared to wild type 
L3 stage worms (Figure 4A), and tbx-9(RNAi) performed 
in a sensitized lin-39 background caused a defect in 
VPC fusion at the L2 (Table 4). tbx-9 RNAi treatment 
did not effect lin-39::GFP expression in the cells P5-P8 
in the embryo (Table 2). In summary, we found that 
TBX-9 may bind multiple sites within a 257 bp fragment 
from lin-39 intron 2, and TBX-9 acts as a positive regu- 
lator of lin-39 expression in the VPCs in the larva. 

The Zinc-finger protein BED-3 binds to site in lin-39 intron 
2 and positively regulates lin-39 expression 

bed-3 encodes a BED zinc-finger protein that is ex- 
pressed in most hypodermal cells, including the seam 
cells and the progeny of the VPCs at the time of L3/L4 
molt [68]. In bed-3 mutants, the granddaughters of P5.p, 
P6.p and P7.p often fail to divide, which suggested that 
BED-3 acts late during vulval induction in the terminal 
divisions of the induced VPCs [68]. In yeast, BED-3 
bound YF10, a 319 base pair fragment from the second 
lin-39 intron that contains ECRs 24-26 (Additional file 1: 
Figure SI and Additional file 2: Figure S2). Bacterially- 
expressed BED-3 protein bound to YF10 in vitro, however, 
we found that BED-3 protein also bound to several other 
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Figure 5 NHR-43, LIN-26, ELT-6 and EGL-18 are necessary for lin-39::GFP expression in the embryo. A - D) GFP expression in embryos 
derived from smg-l; dels4(lin-39::GFP) animals treated for control RNAi (empty vector, panel A) or for RNAi against transcription factor genes 
B) nhr-43; C) lin-26 or D) e/f-6. E - H) GFP expression in embryos carrying pJW3.9::GFP in a wild-type background (E) or in animals carrying 
mutations in e/f-6 (F) or egl-18 (G and H). Embryos shown are at the 'bean' stage of embryogenesis (-360 minutes). All photos were taken under 
identical exposure settings. Note that before the individual P cells interdigitate at the ventral midline, the cells are referred to by their possible 
fates (i.e., P5/6 L and P5/6R), and there are two GFP-expressing cells on the left side and two on the right. 



unrelated DNA fragments, suggesting the purified BED-3 
protein may show non-specific binding in vitro (data not 
shown). 

Despite our inability to validate the BED-3 binding re- 
sult in vitro, we observed that bed-3 RNAi caused a 
strong decrease in the number of animals with wild type 
levels of lin-39::GFP expression in P5.p-P8.p at the L3 
stage (Figures 2G and 3F), and this decrease was also 
seen by qRT-PCR analysis on bed-3(sy702) L3 larvae 
(Figure 4A). Consistent with the decrease in lin-39 ex- 
pression in vivo, bed-3(RNAi) caused a defect in VPC 
fusion at the L2 in a sensitized lin-39 background 
(Table 4). However lin-39::GFP expression did not 



change in embryos derived from mothers treated for 
bedS(RNAi) (Table 2). Therefore, although we could not 
localize a binding site for BED-3 beyond the lin-39 in- 
tron 2 fragment used in the yeast screen, our in vivo 
data indicate that BED-3 is likely to function as a posi- 
tive regulator of lin-39 expression in the vulval precursor 
cells in the larva. 

The GATA factor ELT-6 binds to a lin-39 enhancer that 
directs expression in the P cells in the embryo 

We previously described a 340 bp cis-regulatory element 
from lin-39 that is sufficient to drive GFP expression in 
the embryo in P5-P8, cells which divide to generate lin-39 
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Table 4 Reduction of transcription factor function in a 
lin-39 sensitized background affects VPC fusion 


Strain 


% WT VPCs 


Iin-39(n709ts) FV(RNAi) 


96% 


Iin-39(n709ts) nhr-43(RNAi) 


88% 


Iin-39(n709ts) alr-l(RNAi) 


72%* 


Iin-39(n709ts) ztf-17(RNAi) 


80% 


Iin-39(n709ts) lin-26(RNAi) 


65%* 


Iin-39(n709ts) tbx-9(RNAi) 


72%* 


Iin-39(n709ts) bed-3(RNAi) 


72%* 


Iin-39(n709ts) elt-6(RNAi) 


56%* 



Iin-39{n709ts); wls79(ajm-1::GFP) animals were grown for two generations at 25° 
on f. coli expressing dsRNA targeting the indicated gene (except for bed-3 and 
lin-26 which were grown one generation due to lethality). Larvae at the mid L2 
stage (17 hrs after feeding of newly hatched Lis) were examined for ajm-1:: 
GFP expression in the vulval precursor cells. Wild type animals have either 5 or 
6 VPCS with ajm-T::GFP expression at this time (due to fusion of P3.p with 
hyp7). The percentage of animals showing the wild type pattern is shown (n > 
25 animals). '*' indicates P < 0.05 compared to feeding vector (FV) control 
(Fisher's exact test). 



expressing VPCs and neuroblasts in the ventral midbody 
region (construct p/W3.9, [47]). A 24 bp sequence (SI) in 
the pJW3. 9 enhancer is conserved between C. elegans and 
C. briggsae, and mutation of this site abolished embryonic 
expression from the pJW3.9 reporter [47]. Site SI contains 
the sequence TGATAA, a predicted binding site for a 
GATA family transcription factor, which prefer the motif 
WGATAR [69]. Intriguingly, we found that the transcrip- 
tion factor ELT-6 bound to the p/W3. 9 enhancer fragment 
in our yeast screen (Additional file 1: Figure SI and 
Additional file 2: Figure S2). elt-6 encodes a 367 amino 
acid GATA transcription factor expressed in certain neu- 
rons and hypodermal cells, particularly the seam cells and 
VPCs [55,70]. We performed directed yeast one-hybrid as- 
says with eight other C. elegans GATA factors and found 
that ELT-6 was the only GATA factor that interacted with 
the pJW3.9 enhancer fragment in yeast (Additional file 4: 
Table S5). When the GATA site in pJW3.9 was mutated in 
the yeast reporters, ELT-6 no longer interacted with the 
DNA fragment in yeast (Figure 6A). ELT-6 protein 
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Figure 6 ELT-6 interacts with pJW3.9 through a GATA binding site. A) Yeast strains containing HIS3 and lacZ reporters transformed with 
empty vector (Control) or ELT-6::GAL4AD plasmids (ELT-6) were diluted and replica plated to control (left), 10 mM 3aminotriazole (3AT, middle), 
and XGal (right) plates. Reporters had inserts of either JW3.9 (top) or JW3.9 in which GATA binding site S1 was mutated from TGATAA to GGTACC 
('Mut' bottom). Mutation of the GATA site abolishes the interaction with ELT-6 based on lack of growth on 3AT and lack of blue color on XGal.; 
B) Increasing concentrations of ELT-6 interact with a 40 bp fragment around GATA site S1 (lanes 2-6). Mutation of the GATA site SI from TGATAA 
to GGTACC (M1) abolishes interaction with ELT-6 (lane 8); C) 50 ng of cold wild type site SI (W) can compete effectively for binding of ELT-6 to 
SI (compare lanes 2 and 3). Mutation of GATA site S1 (Ml) reduces but does not abolish the ability to compete (lane 4). Mutation of both site S1 
and a second GATA site on the oligonucleotide (M2) drastically reduces the ability of the oligonucleotide to compete for ELT-6 binding (lane 5). 
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purified from E. coli bound the SI site in vitro, but did not 
bind when the GATA sequence was mutated (Figure 6B). 
ELT-6 binding to SI was abolished when competed with 
excess wild type cold SI probe, but when the SI GATA 
site was mutated, the resulting oligonucleotide (Ml) com- 
peted less well for ELT-6/S1 binding (Figure 6C). When a 
second GATA site at the edge of the 40 bp oligonucleotide 
was also mutated, the ability of the mutated oligonucleo- 
tide (M2) to compete was greatly reduced (Figure 6C). 

To examine regulation of lin-39 expression by ELT-6 
in the embryo, we assayed expression of a reporter con- 
struct containing the wild type 340 bp element, pJW3.9:: 
GFP, in elt-6(gk723) mutant animals. gk723 is an allele 
with 457 bp deletion covering the first intron and sec- 
ond exon of elt-6 and is a presumed null mutation [57]. 
Only 56% of embryos from elt-6(gk723) mutant animals 
displayed a wild type pattern of p/W3.9::GFP expression, 
compared to 94% of control embryos (Figure 5F, Table 3). 
A decrease in the penetrance of expression was also seen 
in the embryos derived from pJW3.9::GFP mothers treated 
with elt-6 RNAi (78%; Table 3). We also examined the ef- 
fect of reduction of elt-6 function on embryonic expres- 
sion of the large lin-39:.GFP reporter that we used to 
assay other transcription factors. When elt-6 RNAi was 
performed on lin-39:.GFP hermaphrodites, only 82% of 
embryos showed the wild type level of GFP expression in 
P5- P8 (Figure 5D, Table 2). Taken together, these results 



indicate that the GATA factor ELT-6 is necessary for 
proper expression of lin-39 in P5-P8 in the embryo, most 
likely via binding to the conserved GATA site in the 
pJW3.9 enhancer, which was previously shown to be ne- 
cessary for enhancer driven GFP expression in the embry- 
onic P cells [47]). 

The GATA factor EGL-18 also regulates lin-39 expression 
in the P cells in the embryo 

The elt-6 open reading frame begins less than 600 bp 
downstream from the end of another GATA factor gene, 
egl-18 (see Figure 7B), and these two genes are tran- 
scribed dicistronically in some tissues [55,70]. The DNA 
binding domains of EGL-18 and ELT-6 are similar, and 
the two genes show genetic redundancy during fate 
specification of the hypodermal seam cells and VPCs 
[55,70,71]. In particular, reduction of function for both 
elt-6 and egl-18 in the larva causes the VPCs to adopt in- 
appropriate cell fates and fuse with the hypodermal syn- 
cytium [55], a phenotype also seen with reduction of 
lin-39 function [28,45]. 

Although EGL-18 did not bind to the 340 bp pJW3.9 
fragment in yeast assays (Additional file 4: Table S5), the 
known functional redundancy of elt-6 and egl-18 and 
their dicistronic expression in hypodermal cells led us to 
test for a role for egl-18 in the regulation of lin-39 ex- 
pression. When the pJW3.9::GFP reporter was moved 
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Figure 7 Direct transcriptional regulators of lin-39 in the embryo and larva. A) Horizontal lines represent 20 kb of genomic DNA 
surrounding the lin-39 locus. The lin-39 transcript is shown below the top line, with boxes representing exons. The next horizontal line shows 
evolutionarily-conserved regions (ECRs; thin vertical lines), the PCR fragments used in the yeast one hybrid assays containing the ECRs (boxes 
labeled 1-12), and two fragments (pJW3.9 shown, JW5 unlabeled) identified previously using an enhancerless GFP assay [47]. Transcription factors 
that bind the lin-39 gene are shown above the line (previously reported) or below the line (reported in this work). B) Model for positive feedback 
loop between egl-18/ elt-6 and lin-39. EGL-18 and ELT-6 act via the GATA site in enhancer pJW3.9 to facilitate initiation of lin-39 expression in the 
embryo, and then LIN-39 acts to positively regulate egl-18/elt-6 expression via the Hox/Pbx binding site in the intron of egl-18 [55]. 
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into two different egl-18 mutant strains, egl-18(ga97) 
and egl-18(nl62) [55], we found that GFP expression in 
P5 - P8 in the embryo decreased from 94% for control 
animals to 60% in egl-18(ga97) and 68% in egl-18(nl62) 
(Figure 5G and H, Table 3). These results suggest that 
EGL-18 is also required for proper lin-39 expression in 
the embryo. To test for redundant function in lin-39 
regulation, we performed RNAi for one gene in the 
background of a mutant for the other (because the genes 
are adjacent, we were unable to build egl-18; elt-6 double 
mutant animals). In neither case was there a significant 
decrease in the number of embryos showing p/W3.9:: 
GFP expression compared to the egl-18 or elt-6 mutant 
strain treated with control vector (Table 3). Thus we did 
not obtain evidence for functional redundancy of elt-6 
and egl-18 in regulating pJW3.9::GFP expression in the 
embryo, even though reduction of function for either 
gene alone affects embryonic GFP expression. 

ELT-6 and EGL-18 also regulate lin-39 expression 
post-embryonically 

We also examined whether elt-6 affects lin-39 expres- 
sion post-embryonically. elt-6 RNAi was performed on 
lin-39::GFP LI larvae and we observed GFP expression 
in the VPCs in the same animals at the L3 stage. A sta- 
tistically significant decrease in the number of VPCs (P5. 
p - P8.p) showing a wild type level of lin-39::GFP expres- 
sion in the VPCs at the L3 stage was seen for elt-6(Ll 
RNAi) animals (Figures 2H and 3G), indicating that elt-6 
also positively regulates lin-39 in the larval VPCs. EGL- 
18 is also likely to regulate lin-39 expression postem- 
bryonically, since qRT-PCR on animals in which both 
egl-18 and elt-6 function is compromised had lower ex- 
pression of lin-39 than elt-6(gk723) mutant animals 
alone (Figure 4B). 

To test the hypothesis that lin-39 is a downstream tar- 
get of ELT-6 and/or EGL-18 in the larva, we asked if 
overexpression of the GATA factors could increase lin- 
39 levels. We overexpressed egl-18 and elt-6 using the 
heat shock promoter and assayed lin-39 levels by qRT- 
PCR for lin-39 one hour after a single heat shock at the 
L2/L3 molt. We found that lin-39 levels went up 1.8 fold 
when ELT-6 was over-expressed, 1.6 fold when EGL-18 
was over-expressed, and 1.7 fold when both ELT-6 and 
EGL-18 were over expressed (Figure 4C), supporting the 
hypothesis that lin-39 is a downstream target of both 
ELT-6 and EGL-18 in the larval VPCs. 

Discussion 

The C. elegans Hox gene lin-39 functions in the mid- 
body region of the developing C. elegans larva, where it 
is expressed in the P cells and their descendants, includ- 
ing the hermaphrodite vulval precursor cells (VPCs) 
[14]. During vulval development, lin-39 expression is 



regulated by Wnt and Ras signaling pathways to facili- 
tate VPC fate specification [26,28,45,46]. To further 
understand the function of lin-39 in hermaphrodite vul- 
val development, we wish to identify cis-acting sites and 
trans-acting factors required for initiation, maintenance 
and regulated expression of lin-39 in the P cells and 
VPCs. Previously, using an enhancerless GFP reporter 
assay, we identified a lin-39 enhancer fragment (JW5) 
that directs expression in the VPC P6.p at the time of 
vulval induction, which responds to Ras pathway activity, 
and which is bound by the Ras pathway effectors LIN-1 
and LIN-31, as well as by LIN-39 itself [47]. We also 
identified a lin-39 enhancer (JW3.9) that directs expres- 
sion in the P cells in the embryo and in their larval de- 
scendants (the VPCs and ventral cord neurons), and 
showed that expression in the embryo depended on an 
evolutionarily conserved site that contains within it a 
putative GATA transcription factor binding site. 

Here we took a complementary approach to identify 
additional factors that may regulate lin-39 expression 
during vulval development. We concentrated on 27 
short, evolutionarily conserved regions (ECRs) that we 
previously identified in the lin-39 gene in C. elegans, 
C. briggsae and C. remanei [47], on the assumption 
that some may represent binding sites for phylogenet- 
ically conserved transcriptional regulators. We used 
the yeast one-hybrid technique to identify transcrip- 
tion factors that could bind to DNA fragments con- 
taining one or more of these ECRs, circumventing the 
requirement that the DNA site be sufficient to direct 
in vivo reporter expression [53]. Having found mul- 
tiple transcription factors that interact with lin-39 
DNA fragments, we looked for effects on endogenous 
lin-39 transcript levels and on expression from a lin- 
39..GFP reporter in vivo when the function of these 
proteins was reduced. Using this approach we identi- 
fied six transcription factors that bind lin-39 promoter 
sequences in yeast and in vitro and regulate lin-39 ex- 
pression and/or function in vivo (Figure 7A, Table 1). 
Three factors, the orphan nuclear receptor NHR-43, 
the hypodermal cell fate regulator LIN-26 and the 
GATA factor ELT-6 positively regulate lin-39 expres- 
sion in the embryo. In addition to those three factors, 
we identified the zinc finger proteins ZTF-17 and 
BED-3 and the T box factor TBX-9 as positive regula- 
tors of lin-39 expression in the larval VPCs. Before 
this work, only four transcription factors were known 
to directly bind at the lin-39 gene and regulate its ex- 
pression; LIN-1, LIN-31, LIN-39 and TRA-1, all of 
which act during larval life (Figure 6B) [37,47]. There- 
fore, by this approach combining phylogenetic conser- 
vation and yeast one-hybrid screening, we have more 
than doubled the number of factors known to directly 
bind to and regulate lin-39, and for the first time we 
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have identified factors that regulate expression of this 
Hox gene in the embryo. 

We believe these factors regulate lin-39 in vivo be- 
cause we observed a reproducible change in expression 
of the full length lin-39::GFP reporter and/or a decrease 
in endogenous lin-39 transcript levels when the function 
of each of these factors was reduced, and in several cases 
we observed a phenotype in a LIN-39 regulated process 
(prevention of VPC fusion). Given their binding to sites 
from the lin-39 genomic region, the simplest model 
would be that these factors directly regulate lin-39 ex- 
pression in the cells we examined (Figure 7A). However, 
in the current work we did not pursue in vivo binding 
studies for any of these factors; although published re- 
sults from the modENCODE project show binding of 
ALR-1 to the evolutionarily conserved site we identified 
at the relevant time in development [40,63]. Addition- 
ally, while existing GFP reporters for four of these pro- 
teins (TBX-9, NHR-43, ALR-1 and ZTF-17) are turned 
on in some &-39-expressing cells, expression in the P cells 
or VPCs has not been directly observed [57,59,61,64,67]. 
GFP reporters suffer from the caveat that all the elements 
required to recapitulate endogenous expression may not 
be present, or the expression may be weak or dynamic. 
However, until we can verify expression of each transcrip- 
tion factor in the embryonic P cells or larval VPCs, and 
show evidence of binding to the sites we identified in vivo, 
it remains possible that some of these factors regulate lin- 
39 indirectly via another transcription regulator, or even 
act non cell-autonomously on lin-39 expression. The 
strongest case for direct regulation can be made for BED- 
3, which regulates lin-39 expression during larval life, and 
for LIN-26 and ELT-6, which regulate lin-39 expression in 
the embryonic P cells. 

BED-3 

bed-3 encodes a 599 amino acid zinc finger protein with 
a BED DNA binding domain, and bed-3 mutants show 
an Egl- laying (Egl) defective phenotype and defects in 
the terminal cell divisions of the descendants of the 
VPCs [68]. An intronic enhancer element from bed-3 di- 
rects GFP expression in VPC descendants, leading to the 
hypothesis that BED-3 functions in the terminal cell di- 
visions of vulval cells during vulval induction [68]. We 
found that bed-3(RNAi) showed the strongest effect on 
full length lin-39::GFP expression in the VPCs; in wild 
type animals the VPCs P5.p - P8.p showed lin-39::GFP 
expression 89% of the time (averaged over all four cells), 
while in bed-3(RNAi) animals these cells showed GFP 
expression only 37% of the time. This suggests that 
BED-3 could be acting at an earlier stage in vulval devel- 
opment. Consistent with this, we observed that the bed- 
3 enhancerr.GFP reporter does show expression in the 
VPCs before they divide (Additional file 8: Figure S7), 



which would be consistent with BED-3 acting upstream 
of lin-39 in the VPCs themselves. Therefore, although 
we could not verify BED-3 binding to the YF10 site 
in vitro with purified protein, we believe BED-3 is likely 
to be a positive regulator of lin-39 in the larval VPCs be- 
fore vulval induction, as well as functioning in the subse- 
quent cell division of their progeny, as previously 
reported [68]. 

LIN-26 

LIN-26 is a 438 amino acid zinc finger transcription fac- 
tor that shows continuous expression in all embryonic 
and larval hypodermal cells after their birth [65]. Loss of 
lin-26 causes hypodermal cells to adopt incorrect cell 
fates or degenerate after their birth, resulting in embry- 
onic lethality [65,66]. Ectopic expression of LIN-26 in 
the early embryo can induce cells to adopt hypodermal- 
like fates [72]. These results suggest that LIN-26 is re- 
quired to specify and/or maintain the hypodermal cell 
fate [65] . lin-26 itself is positively regulated by the GATA 
factor ELT-1, which is another global regulator of hypo- 
dermal cell fate [73]. We found that lin-26 RNAi caused 
a weak but significant decrease in the expression of the 
full length lin-39::GFP reporter in the P cells in the em- 
bryo; 100% of wild type embryos show GFP expression, 
compared to 86% of lin-26(RNAi) animals. LIN-26 binds 
in vitro to a 51 bp sequence located approximately 2 kb 
upstream of the lin-39 start codon (Figure 6B). Recent 
data recording reporter gene expression from live devel- 
oping embryos shows that a Un-26::mCherry transcrip- 
tional reporter shows expression in the mothers of the P 
cells shortly before they divide, at the same time as these 
cells also begin to show expression of a lin-39::mCherry 
reporter [74] (Additional file 9: Figure S8). Based on 
these results, and the known function of LIN-26 in hy- 
podermal cell fate, we propose that LIN-26 positively 
regulates lin-39 expression in the embryonic P cells, and 
that the function of this regulation is to aid in initial P 
cell fate specification and/or to maintain the P cell iden- 
tity once established. 

ELT-6 

elt-6 encodes a 367 amino acid GATA family transcrip- 
tion factor. The elt-6 gene is immediately downstream 
from another GATA factor gene, egl- 18 (Figure 7B), and 
reporter gene experiments suggest elt-6 may be expressed 
on its own and as part of a dicistronic message with egl- 18 
[70]. These genes are expressed in many cells in the em- 
bryo, including the descendants of the MS and AB blasto- 
meres (which give rise to the P cells). In the larva, they are 
expressed strongly in the lateral hypodermal seam cells 
and weakly in the VPCs [55,70]. These two GATA factors 
share 76% identity in their DNA-binding domains, and 
have been shown to act redundantly in seam cell 
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development in both the embryo and larva [70,71]. We 
previously showed that a 340 bp lin-39 enhancer (JW3.9) 
directs GFP expression in the embryonic P cells, and a 
conserved GATA site in the enhancer was necessary for 
expression [47] . Here we identify the GATA factor ELT-6 
as binding to this lin-39 enhancer in yeast and in vitro, 
and show that binding was dependent on the GATA se- 
quence. elt-6 RNAi showed a decrease in expression of 
the full length lin-39::GFP reporter in embryonic P cells, 
and this result was recapitulated with the pJW3.9 reporter 
in elt-6(gk723) mutant and elt-6(RNAi) embryos. These 
data suggest that ELT-6 is required for proper expression 
of the Hox gene lin-39 in the embryonic P cells. Consist- 
ent with this, recent data from live recordings of develop- 
ing embryos shows that elt-6 expression in the mothers of 
P3 - P8 begins before lin-39 expression is first seen ([74]; 
Additional file 8: Figure S7). This is the first example to 
our knowledge of a phenotype caused by reduction of elt- 
6 function alone, elt-6 RNAi treatment of newly hatched 
LI larvae also led to a weak reduction in lin-39::GFP ex- 
pression in the L3 VPCs, and overexpression of ELT-6 or 
EGL-18 increased endogenous lin-39 expression in larvae. 
These results suggest that ELT-6 and/or EGL-18 may con- 
tinue to positively regulate lin-39 expression in larval life 
in the VPCs. 

Given these results, it is interesting that we previously 
showed, in collaboration with the Rothmann laboratory, 
that the egl-18/elt-6 locus was likely to be a downstream 
target of LIN-39 in the larval VPCs [55]. In that work an 
800 bp enhancer element was identified in intron 2 of 
the egl-18 gene that directs GFP expression in the VPCs 
and their descendants starting in the L2 stage. This en- 
hancer contains two Hox protein-binding sites that are 
bound in vitro by LIN-39 and its binding partner CEH- 
20, and mutation of one site eliminated enhancer-driven 
GFP expression. Finally, overexpression of egl-18 from 
the heat shock promoter was able to partially rescue vul- 
val defects in lin-39(RNAi) animals. These data led to 
the model that egl-18lelt-6 is a downstream target of 
LIN-39 during vulval development. 

Combining these previous data with our current re- 
sults, our working model suggests LIN-26 and ELT-6 are 
involved in the initiation of expression of the Hox gene 
lin-39 in the P cells in the embryo (Figure 7A). lin-39 
may also be regulated by EGL-18 at this time, since 
egl-18 mutants had reduced enhancer-driven GFP ex- 
pression in the P cells, and EGL-18 protein bound to the 
p/W3.9 GATA site in vitro (W. Liu and D. Eisenmann, 
unpublished results). Once the fate of these cells is 
established, a positive feedback loop is established be- 
tween ELT-6, acting via the upstream lin-39 enhancer 
pJW3.9, and LIN-39, acting via the intronic enhancer in 
egl-18 [55]. We propose that this positive feedback loop 
helps maintain expression of these genes and the fate of 



these cells and their descendants, the VPCs, during sub- 
sequent embryonic and larval development (Figure 7B). 
Interestingly, we have seen another example of lin-39 
feedback regulation. The zinc finger protein SEM-4 was 
previously shown to positively regulate lin-39 expression 
in VPCs [49]. We have found that in animals overex- 
pressing either LIN-39 or LIN-39 and CEH-20, sem-4 
expression is decreased, suggesting that feedback mecha- 
nisms to decrease lin-39 levels when they are elevated 
may also exist (J. Siegel and D. Eisenmann, unpublished 
results). Recent chromatin immunoprecipitation experi- 
ments from larvae verify binding of tagged LIN-39 pro- 
tein upstream and in introns of both the egl-18 and sem- 
4 genes [40,63]. 

Conclusion 

While much is known about the initiation and regula- 
tion of Hox gene expression in Drosophila and verte- 
brates, less is known outside of these well-studied 
groups. Our laboratory and others have been studying 
the expression and function of the Hox gene lin-39 dur- 
ing vulval development in the nematode C. elegans. We 
used the yeast one-hybrid (Y1H) method to identify 16 
transcription factors that interact with specific regions of 
the lin-39 gene, and further characterized several factors 
(ALR-1, BED-3, ELT-6, LIN-26, NHR-43, TBX-9 and 
ZTF-17) showing that their function is required for 
proper expression of lin-39::GFP reporters and endogen- 
ous lin-39 in vivo. This work greatly expands the number 
of factors known to directly regulate lin-39 expression. 
Given the known caveats of the yeast one-hybrid tech- 
nique (absence of specific posttranslational modifications, 
lack of heterodimeric binding factors), and our emphasis 
on characterization of factors with known expression or 
phenotypes in the vulval cells, it is likely that there are 
additional factors that regulate lin-39 in the embryo and 
larva that we did not identify in our screens. Given the im- 
portant role of Hox genes in patterning the developing 
metazoan body, it is not surprising that Hox gene expres- 
sion is found to be complicated in those species where it 
has been examined closely. Our results suggest that ex- 
pression of the C. elegans Hox gene lin-39 in the P cells 
and VPCs during vulval development may be regulated by 
a large number of transcription factors, each making a 
small contribution to overall lin-39 expression on its own. 
This is consistent with in vivo data from C. elegans show- 
ing that the average worm gene is bound by several tran- 
scription factors at one time in larval development 
[40,63]. Such a mechanism may ensure a robustness of ex- 
pression for important developmental regulators like Hox 
genes. Transcription factors behaving in this manner 
would also not be identified in genetic screens, since each 
makes a small overall contribution to lin-39 expression, 
and feedback mechanisms may exist to compensate for 
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reductions in lin-39 transcript levels. Finally, for the first 
time, we have identified factors required for lin-39 expres- 
sion in the embryo (NHR-43, LIN-26 and ELT-6), and our 
results with EGL-18 and ELT-6, combined with our earlier 
work on LIN-39 regulation of egl-18lelt-6, hint at a posi- 
tive feedback mechanism to maintain lin-39 expression in 
the vulval lineages. The identification of LIN-26, ELT-6, 
and NHR-43 will help us further characterize the mecha- 
nisms for the initiation of Hox gene expression in the 
nematode, allowing us to make comparisons across meta- 
zoan phyla about the mechanisms utilized to regulate this 
essential class of developmental regulators. 

Methods 

Yeast one-hybrid (Y1 H) assays 

Y1H assays were performed using two methodologies: a 
traditional library 'transformation' screen, and a robotically- 
assisted 'mating' screen ("enhanced" Y1H, eYIH). Briefly, 
each DNA fragment of interest (YF1 - YF12, p]W3.9; see 
Additional file 1: Figure SI) was cloned into reporter vec- 
tors pMW#2 (HIS3) and pMW#3 (LacZ) and the resulting 
constructs were sequentially integrated into the genome 
of yeast strains BY5444 and YM4271 to generate thirteen 
"DNA bait" strains. BY5444 and YM4271 are isogenic for 
multiple marker genes (MATa um3-52 his3-200 ade2-101 
lys2-801 leu2-3,112 trpl-901 tyrl-501 gal4-A512 gal80- 
A538 adeSr.hisG) but BY5444 does not mate efficiently 
with the Yalphal867 strain used in eYIH assays. 

The 'transformation' screen was performed as de- 
scribed [53]. Each of the thirteen bait strains (BY5444 
background) was transformed with DNA from a 
commercially-available, C. elegans transcription factor li- 
brary containing 755 plasmids that each express one C. 
elegans transcription factor fused in frame with the yeast 
GAL4 activation domain (Open Biosystems; [75]). Po- 
tential interactions were identified as those colonies that 
turned blue on plates containing X-Gal, and which grew 
on plates containing a higher concentration of 3- 
aminotriazole than the control strain grew on. The 
eYIH screens were performed as described [54]. Each of 
the thirteen bait strains (both BY5444 and YM4271 
backgrounds) was mated with a collection of 936 YlHal- 
phal867 strains of the opposite mating type, each of 
which expresses a single C. elegans transcription factor 
fused to the GAL4 activation domain. Matings were 
done in quadruplicate. Potential interactions were identi- 
fied as those for which at least two of the four colonies 
exhibited higher expression of both reporters than con- 
trol yeast, as assayed by blue color and growth on SC- 
His + 5 mM 3-aminotriazole + X-Gal. For both methods, 
plasmids were recovered from positive yeast, sequenced 
to identify the C. elegans gene, and then retransformed 
back into the appropriate yeast strain to confirm the 
interaction. From both screens, only those interactions 



that repeated after retransformation were considered 
true positives. For the haploid transformation screen, 
while most baits gave many 3-aminotriazole resistant 
colonies, only RF1, RF2, RF4 and RF6 gave double posi- 
tive colonies, and only interactions with NHR-43, ALR-1 
and ZTF-17 repeated (all three factors were identified 
multiple times). For the robotically assisted mating 
screens, the results are presented in Additional file 4: 
Table S2-S4. 



Genetic methods, alleles and strains 

Methods for culture and genetic manipulation of C. ele- 
gans were performed as described [76]. Bristol variety 
(strain N2) of C. elegans was used as wild-type. Experi- 
ments were performed at 20°C unless noted. Genes and 
alleles used in this work are described in [77] and 
Wormbase [57]. 
LGII: rrf-3(pkl426), Un-26(nS6) 

LGIII: tbx-9(ok2473), tbx-9(ms31), unc-119(ed3), pha-1 
(e2123), lin-39(n709ts) 

LGIV: nhr-43(tml381), bed-3(sy702), bed-3(sy705), elt- 
6(gk723), egl-18(ga97), egl-18(nl62), egl-18(ok290) 
LGX: alr-l(ok545) 
Strains used: 

lin-39::GFP: smg-l(el228); him-5(el490); dels4[lin- 

39TN::GFP; dpy-20(+); ajm-l::GFP] [46] 

nhr-43::gfp: unc-119(ed3); Ex[C29E6.5::gfp; unc-119 (+)] 

[64] 

bed-3::GFP: unc-119(ed3); syEx962[pTI06.29; unc-119 
(+); pBSKSII(+)] [68] 

pJW3.9::GFP: pha-l(e2123); deExl05[pJW3.9::gjp; pha- 
1(+); ajm-l::GFP] [47] 
Strains created for this work: 

hs:xontrol: unc-119(ed3); deExl06[pPD49.78; unc-119 
(+); ajm-l::GFP] 

hs::elt-6: unc-119(ed3); deExl07[pLGl; unc-119(+), 
ajm-h-.GFP] 

hs::egl-18: unc-119(ed3); deExl08[pPK8; pPK9; unc-119 
(+); ajm-h-.GFP] 

hs::egl-18+ hs::elt-6: unc-119(ed3); deExl09 [pLGOl; 
pKK8;pKK9 unc-119 (+); ajm-h-.GFP] 

RNA interference 

RNA interference (RNAi) was performed using the bac- 
terial 'feeding' method in which dsRNA for the gene of 
interest is produced in E. coli strain HT115 and ingested 
by worms [78]. For feeding of LI larvae, eggs from 
strains to be tested were placed on plates without E. coli 
for at least 18 hours at 20°C, then semi-synchronized LI 
larvae were washed off, placed onto plates with the de- 
sired HT115 RNAi strain and grown at 20°C to the ap- 
propriate stage for scoring. For feeding of PO animals, LI 
stage PO larvae were put onto the desired RNAi plates, 
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grown to adulthood and their Fl progeny were analyzed 
at the appropriate stage. 

Reporter gene analysis 

Worms carrying GFP reporter constructs were analyzed 
using fluorescence microscopy on a Zeiss Axioplan 2 at 
the desired developmental time. GFP expression in live 
animals was captured using a Nikon DXM 1200 digital 
camera and the ACT-1 program (version 2.12). For dels4 
[lin-39::GFP] and pJW3.9::GFP reporter analysis, the per- 
centage of animals showing an intensity of expression 
similar to wild-type was determined. For embryonic ex- 
pression in the cells P5 - P8 the number of embryos 
showing wild type expression in all four cells was re- 
corded. For larval expression in the cells P5.p - P8.p the 
intensity of GFP expression in the individual cells at the 
L3 stage was analyzed using ImageJ [79] and pixel 
counts for each cell were recorded (after subtraction of 
background). For each RNAi treatment, at least 20 ani- 
mals at the L3 stage were photographed, under identical 
conditions. Expression data was statistically analyzed 
using an unpaired i-test and both P values and SD were 
gathered. All dels4 experiments were carried out in a 
smg-1 mutant background, in which nonsense-mediated 
decay is abrogated, leading to more robust reporter ex- 
pression [46]. 

Creation of strains expressing EGL-18 and ELT-6 from the 
heat shock promoter 

For heat-shock induced expression of EGL-18, we used 
the previously described constructs pKK8 and pKK9 in 
which the entire coding region of egl-18 is inserted into 
the heat shock promoter vectors pPD49.78 and pPD49.83 
respectively [55]. A similar heat shock expression con- 
struct, pLGl (gift of L. Gorrepati), was made by cloning 
elt-6 genomic DNA from the start to stop codon into vec- 
tor pPD49.78 (available at [80]; gift of A. Fire, Stanford 
University School of Medicine, Stanford, CA). To create 
transgenic strains, DNA for these constructs was co- 
injected at 25 ng/ul with unc-119 (+) (100 ng/ul) and 
ajm-l::gfp DNA (50 ng/ul) [81] into unc-119(ed3) her- 
maphrodites. For the hsr.control strain, the empty heat-shock 
vector pPD49.78 was used. For the hs::elt-6 + hs::egl-18 
strain, pLGl, pKK8 and pKK9 were injected together. 
Several lines rescued for the Unc phenotype were recov- 
ered for each injection and the ones with the highest 
percentage transmission to progeny were analysed. 

Heat shock protocol 

Embryos from transgenic animals containing the arrays 
hsr.control, hs::elt-6, hs::egl-18, and hs::elt-6 + hs::egl-18 
were collected from hypochlorite-treated gravid animals 
and hatched on NGM plates without food for at least 
18 hours at 20°C, allowing for early LI stage arrest. 



Synchronized LI animals were fed with OP50, grown to 
the early L3 stage (25 hours post feeding at 20°C), heat 
shocked at 33°C for 30 minutes, then transferred back to 
20°C and collected 30 minutes later for RNA isolation. 

RNA isolation 

For genes with existing mutant strains, animals were 
synchronized and grown to early L3 stage (26 hours after 
feeding at 20°C) to be collected. For genes without exist- 
ing mutant strains, RNAi was performed on wild type 
worms. For heat shocked worms, animals were treated 
with the heat shock program described above. Three bio- 
logical replicates were performed for each strain and RNA 
was extracted using the RNAeasy mini kit (Qiagen). 

qRT-PCR 

The mRNA fraction of extracted total RNA pools was 
reverse transcribed to cDNA using iScript™ cDNA syn- 
thesis kit (Bio-Rad) and used in triplicate qRT-PCR reac- 
tions run on an iCycleriQ real-time PCR machine 
(BioRad). Relative expression ratios were calculated from 
observed Ct values using the AACt method [fold change = 

tT2 \ targetACt (control-sample)/ (tz \ RefACt (control-sample) rnoll 
^target) / (trefj l°^JJ 

with the house keeping gene gpd-2 as a reference [83]. For 
mutant strains, N2 animals were the control. For RNAi- 
treated worms, N2 animals treated with the empty RNAi 
'feeding vector' (L4440) (available at [80]; gift of A. Fire, 
Stanford University School of Medicine, Stanford, CA) 
were the control. For heat shocked strains, hsr.control ani- 
mals which underwent the heat shock treatment were the 
control. Three biological replicates were used for each 
sample. The data was statistically analyzed using unpaired 
t-test and both P values and SD were gathered with 
Graphpad software. Primers for gpd-2 and lin-39 were 
gpd-2¥W (CCTCTGGAGCCGACTATGTC), gpd-2RV 
(TGGCATGATCGTACTTCTCG), lin-39¥W (CGGAGA 
TCAGTCACTATGCT) and lin-39RV (CCGCGTGAAC 
CTCCTGTAGT). 

Site direct mutagenesis 

Site directed mutagenesis was performed using the 
Quick Change site directed technique following manu- 
facturer's instructions with the high fidelity DNA poly- 
merase Pfu Turbo (Stratagene). Plasmid DNA was 
extracted and sequenced to confirm the presence of 
mutations. 

Protein purification 

Full length cDNAs for nhr-43, alr-1, ztf-17, lin-26, bed- 
3, tbx-9, and elt-6 were obtained by PCR amplification 
from plasmids obtained from the Open Biosystems Y1H 
library [75] and cloned individually into the plasmid 
pQE-80 L (Qiagen), which introduces a 6His tag at the 
N terminus. Plasmids were transformed into E. coli 
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strain BL21, grown to an OD of 0.7 and induced with 
1 mM IPTG for 4 hours. For purification of NHR-43, 
ALR-1, LIN-26, BED-3, TBX-9, and ELT-6, cells were 
centrifuged, resuspended in lysis buffer (50 mM, 
NaH 2 P0 4 , 300 mM NaCl, 10 mM imidazole, pH 8.0 plus 
protease inhibitors [Sigma]), sonicated for 3 minutes on 
a Branson Sonifier 450, and centrifuged. The super- 
natant were collected and loaded on a Ni 2+ -NTA agar- 
ose column (Qiagen), washed five times with wash 
buffer (50 mM NaH 2 P0 4 , 300 mM NaCl, 20 mM imid- 
azole, pH 8.0) and eluted with elution buffer (50 mM 
NaH 2 P0 4 , 300 mM NaCl, 250 mM imidazole, pH 8.0). 
For purification of ZTF-17, cells were centrifuged, resus- 
pended in sonication buffer (20 mM Tris 7.9, 500 mM 
NaCl), sonicated and centrifuged. The pellet was resus- 
pended in solubilizing buffer (6 M Guanidine, 20 mM 
Tris 7.9, 300 mM NaCl, 5 mM Imidazole, 5% glycerol 
plus protease inhibitors [Sigma]). Insoluble material was 
pelleted by centrifugation and the supernatant was 
loaded onto a Ni 2+ -NTA agarose column, washed five 
times with buffer and eluted with elution buffer 2(6 M 
Guanidine, 300 mM NaCl, 300 mM Imidazole, 40% gly- 
cerol, 20 mM HEPES pH7.5). Protein was diluted to 
0.5 mg/ml with elution buffer, dialyzed at 4°C overnight 
against dialysis buffer (20 mM HEPES pH7.5, 300 mM 
NaCl, 40% glycerol) and concentrated to 500 ng/ul using 
a Centricon-10. 
The primers for cDNA amplification were: 
NHR-43: FW:GCGCGGATCCATTAGCGGCCCATTT 
CTTCAC/RV: GCGCGCAAGCTTTTAGATTGAGTA 
CAAGTAGGC 

ALR-1: FW: GCGCGCATGCCCCGAGTTGAAGAAA 
GAAGA/RV: GCGCGCAAGCTTTCATGAACTTTCT 
TCTTTTG 

ZTF-17: FW: GCGCGCATGCCTGCGCTACCAGGC 
GTCCGTG/RV:CGCCCCGGGCTATTTTACTCTAAG 
AAATA 

LIN-26: FW: GCGCGCATGCCTTTCTAAATTTGTG 
GTAGTC/RV: CGCCCCGGGCTACACCAATGGTTG 
AGCCAT 

TBX-9: FW: GCGCGAGCTCTCCAAAGTCAAAGTA 
TCA/RV: CGCCCCGGGTCAACCAACAATATCAAT 
BED-3: FW: GCGCGCATGCCAGACCCAAAGTCCA 
TTT/RV: CGCCCCGGGTCAAACAAGTTGATCAAT 
ELT-6: FW: G CG CG C ATG C ACGTCGTCG AAGG AG 
GAGATG/RV: CGCCCCGGGTCAGGGAGACTTGCG 
CTGCTC 

Electrophoretic mobility shift assays 

For DNA probes larger than 100 bp, DNA fragments 
were produced by PCR, and 5 pmole was labeled with 
32 P using T4 polynucleotide kinase (NEB). For DNA 
probes less than 100 bp, one oligonucleotide was labeled 
with 32 P, then annealed with the complementary 



oligonucleotide. The labeled double stranded DNA 
probes were purified using Centri-Spin-20 columns 
(Princeton Separation), and their activity was measured 
using a standard scintillation counter. The oligonucleo- 
tides used to make EMSA probes are listed in Additional 
file 5: Figure S6. 

DNA binding reaction were set up at 4°C in a volume 
of 20 ul in final buffer conditions of 50 mM KC1, 20 mM 
HEPES, pH 7.9, 0.2 mM EDTA, 0.5 mM DTT, 3 mM 
MgCl 2 , 1 fig poly(dl-dC), 0.5 ug/ul BSA, 5000 cpm 32 P- 
labeled DNA probe and 0-800 ng of protein. Reactions 
were incubated 20 minutes on ice, then loaded onto 
4.5% (probe size smaller than 100 bp) or 6% (probe size 
larger than 100 bp) polyacrylamide gels. Gels were run 
in 0.5xTBE buffer at 150 V for 2.5 hours (4.5% gel), or 
220 V for 2 hours (6% gel), dried and analyzed using a 
Storm 860 phosphorimager (Molecular Dynamics). For 
competitive binding experiments, the indicated concen- 
tration of cold competitor was added to the reaction be- 
fore the protein. 

VPC fusion assay in lin-39 sensitized background 

For nhr-43, alr-1, ztf-17, tbx-9, elt-6 and control RNAi, 
P0 Un-39(n709ts); wk79(ajm-l::GFP) larvae were grown 
to adulthood at 25°C on the desired RNAi plates, then 
semi-synchronized Fl progeny were grown to the late L2 
stage, when ajm-l::GFP expression in the vulval precursor 
cells was examined (17 hours post feeding). For bed-3 and 
lin-26 RNAi, lin-39(n709ts); wk79(ajm-l::GFP) eggs were 
hatched in M9 overnight at 20°C, then semi-synchronized 
LI larvae were washed off, placed onto plates with the 
desired HT115 RNAi strain and grown at 25°C to the late 
L2 stage and scored. At least 25 animals were scored for 
each RNAi treatment using fluorescence microscopy on a 
Zeiss Axioplan 2. Wild type animals have either 5 or 6 
VPCS with ajm-l::GFP expression at this time (due to 
fusion of P3.p with hyp7). The percentage of animals lack- 
ing ajm-1:: GFP expression in P4.p - P8.p was determined. 
The data was statistically analyzed using Fisher's exact test. 
2X2 contingency tables were analyzed using Fisher's exact 
test and P values were gathered. 

Additional files 

* ^ 

Additional file 1: Figure SI. lin-39 genomic region and fragments 
used in yeast one-hybrid screens. The top line shows 26 kb around the 
lin-39 locus with base locations on chromosome II shown. Two lin-39 
transcripts and the upstream microRNA gene mir-231 are diagrammed 
below. JW1-10 (bottom) are lin-39 genomic regions previously used in 
reporter gene analysis; those in green drove GFP expression in vivo in 
lin-39 expressing cells [47], pJW3.9 (orange box) is a 340 bp subfragment 
that drove GFP expression in P5 - P8 in the embryo [47]. Thirty-one 
evolutionarily-conserved regions (ECRs) with >75% identity in the lin-39 
gene from three Caenorhabditis species were previously identified [47]. 
ECRs are shown on the middle line as unlabeled, vertical black lines. 
Twenty-seven ECRS were grouped into 1 2 PCR fragments (numbered 1 -1 2; 
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red boxes). The twelve fragments (YF1-12) and pJW3.9 were used as 'baits' 
in Y1 H screens. 

Additional file 2: Figure S2. Fourteen transcription factors interact with 
lin-39 genomic fragments in yeast one-hybrid assays (strain BY5444). 
Thirteen strains in background BY5444 containing fragments (YF1 - 12, 
pJW3.9) from the lin-39 genomic region were screened for interactions 
with C elegans transcription factors by both the haploid library transformation 
method, and the robotically-assisted mating method (see Methods). Any 
positive interactions from primary screens were retested by retransformation 
of the rescued interacting plasmid back into the appropriate 'bait strain (see 
Additional file 4: Table S2). Shown here are three-fold serial dilutions of each 
strain grown on SC-His-Ura-Trp plates with no 3AT (control plates), SC-His- 
Ura-Trp with 1 0 mM 3AT, and plates with X-gal. Bait strains transformed with 
the pDEST-AD empty vector were the 'Control'. Positive interactions were 
considered those showing more growth on 3AT and/or more blue color on 
XGal than the control. Note that the YF8 strain shows considerable self- 
activation on 3AT, but positive interactions for NHR-1 1 1 and FLH-1 are 
observable on the XGal plates. 

Additional file 3: Figure S3. Four transcription factors interact with 
lin-39 genomic fragments in yeast one-hybrid assays (strain YM4271). 
Thirteen strains in background YM4271 containing fragments (YF1 - 12, 
pJW3.9) from the lin-39 genomic region were screened for interactions 
with C. elegans transcription factors by the robotically-assisted mating 
method (see Methods). Positive interactions from the primary screen 
were retested by retransformation of the rescued interacting plasmid 
back into the appropriate 'bait' strain (see Additional file 4: Table 53). 
Shown here are three-fold serial dilutions of each strain grown on SC- 
His-Ura-Trp plates with no 3AT (control plates), SC-His-Ura-Trp with 
10 mM 3AT, and plates with X-gal. Bait strains transformed with the 
pDEST-AD empty vector were the Control'. Positive interactions were 
considered those showing more growth on 3AT and/or more blue color 
on XGal than the control. 

Additional file 4: Table SI. Sequences of 12 YFs and pJW3.9. Table S2. 
Yeast one-hybrid screen with strain BY5444 (screen YM2). Table S3. Yeast 
one-hybrid screen with strain YM4271 (screen YM2). Table S4. Summary 
of yeast one-hybrid screen results. Table S5. Only GATA factor ELT-6 
binds pJW3.9 in yeast, elegans GATA factors. Table S6. Oligonucleotides 
for EMSA probes. 

Additional file 5: Figure S6. TBX-9 binds a 79 bp region in fragment 
YF9. A) The top line shows fragment YF9 and the location of ECRs 21 - 
23. Subfragments A - D are diagrammed below. Shading indicates the 
fragment bound by TBX-9 in vitro; B) Gel mobility shift assay with TBX-9 
protein purified from £ coli and labeled fragment YF9 (lane 2) and 
competition with subfragments A - D (lanes 3-6). Arrowhead indicates 
free probe; arrow indicates the protein DNA complex that can be com- 
peted by fragment C. Two other bands appear with added TBX-9 protein: 
the upper band is competed away by all four competing fragments, 
while the lower band is competed by none of the fragments. The nature 
of these complexes is unknown, although they are likely to represent 
non-specific binding by TBX-9 or another protein. C) Gel mobility shift 
assay with TBX-9 protein and labeled fragment subfragments A - D. The 
bracket indicates the migration locations of free probes (which differ in 
size, see panel A). Arrow indicates a complex with TBX-9 and subfragment C. 

Additional file 6: Figure S4. ZTF-17 binds to fragment YF4-4 in yeast. 
BY5444 'bait' strains containing subfragments of YF1 (YF4-1 to YF4-4; see 
Figure 1) were made and transformed with the plasmid encoding 
ZTF-1 7::Gal4AD or with the pDEST-AD empty vector as control. Three-fold 
serial dilutions of each strain grown on SC-His-Ura-Trp plates with no 
3AT (control plates), SC-His-Ura-Trp with 10 mM 3 AT, and plates with 
X-gal are shown. Fragment YF4-1 shows strong self-activation on 3AT. 
Only fragment YF4-4 shows a strong interaction with ZTF-17 on both 
3AT and XGal plates. 

Additional file 7: Figure S5. LIN-26 binds to fragment YF4-3 using 
yeast. BY5444 'bait' strains containing subfragments of YF1 (YF4-1 to 
YF4-4; see Figure 1) were made and transformed with the plasmid 
encoding LIN-26::Gal4AD or with the pDEST-AD empty vector as control. 
Three-fold serial dilutions of each strain grown on SC-His-Ura-Trp plates 
with no 3 AT (control plates), SC-His-Ura-Trp with 10 mM 3 AT, and plates 



with X-gal are shown. Fragment YF4-1 shows strong self-activation on 
3AT. Only fragment YF4-3 shows a strong interaction with LIN-26 on 
both 3 AT and XGal plates. 

Additional file 8: Figure S7. bed-3::GFP expression in the larval VPCs. 
Expression from strains carrying syEx962 which contains the bed-3::GFP 
reporter pTI06.29 [68]. GFP expression is seen in the Pn.p cells in the LI 
stage (A) and L2 stage (B), including the VPCs P3.p - P8.p (indicated by 
white bars). 

Additional file 9: Figure S8. Expression of transcription factor genes 
elt-i, elt-3, lin-26 and lin-39 in the ABpra lineage in developing embryos. 
The images shown are taken directly from the Expression Patterns in 
Caenorhabditis web site [84]. As described [74], live images were recorded 
from developing embryos that expressed a histone:mCherry fusion protein 
driven from the upstream promoter sequences for each of the indicated 
genes. Expression levels in individual embryonic cells were characterized 
and diagrammed on the cell lineage chart. An expression scale for all 
experiments is shown at the top and shows the fluorescence intensity of 
the reporter construct in the individual cells; time (in minutes) for each 
experiment is shown along the left. The data shown are from the following 
experiments: 20080128_elt-l_3.html, 2007081 7_e/f-6_5.html, 20080805_//n- 
26_5_Ll.html and 20071 01 5Jin-39_9. 
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